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f"^ ] Abstract. Evidence is presented that the damped Lya absorption sys- 

CN| ' terns are the high-redshift (z > 3) progenitors of galaxy disks. I discuss 

O , , kinematic evidence that the damped Lya systems are rotating disks. I 

QJ ■ also discuss implications of the lack of metal-poor damped Lya systems 

C/^ , with line width Av > 100 km s^^. I then present new evidence stemming 

QQ ' from correlations between element-abundance ratios and [Fe/H], which 

connects damped systems to the thick stellar disk of the Galaxy. I discuss 
the connections between damped Lya systems and Lyman break galaxies, 
and how [CII] 158 fiui emission from damped Lya systems discriminates 
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^vq I among competing theories of galaxy formation. 
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^^ I 1. Introduction 

Q . When did the disks of ordinary galaxies form? This is the central question to be 

in ' addressed in my talk. At this meeting, evidence was presented suggesting the 

t^ . bulk of the disk population was in place hy z ^ 1. While mergers between disks 

and infall of gas onto disks undoubtedly occur subsequently, the surveys by Ellis, 
Lilly, and Faber all rule out substantial changes in the stellar content of disks 
since z = 1. Further evidence against significant changes in the disk populations 
comes from the Galaxy. From her studies of stars in the thick disk, Wyse (this 
a I volume) concluded that the Galaxy has not undergone a major merger during 

the past 12 Gyr. Provided the Galaxy is typical, this adds strong support to 
the idea that disks, in this case the thick disk, were in place before z = 1. 

How long before? This question is crucial since hierarchical cosmologies 
predict that only a small fraction of current disk galaxies were in place by z ~ 
3. Rather, CDM models predict the bulk of the protogalactic mass distribution 
to be in low-mass subgalactic progenitors at these redshifts. This is a generic 
property of CDM cosmogonies, and one that should be tested. 

In my talk I will describe observations that not only test hierarchical the- 
ories of galaxy formation, but also probe the nature of the protogalactic mass 
distribution. The work centers on the damped Lya systems, because of com- 
pelling evidence that these H I layers supply the baryons comprising visible 
stars in current galaxies (cf. Storrie-Lombardi & Wolfe 2000). Therefore, my 
talk focuses on gas rather than stars (but see § 4) . Gas provides a more accurate 
picture of the baryonic content of protogalaxies because most of the baryons at 
high z were gas. 
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2. Kinematics: The Case for Rotation 
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Figure 1. "Edge- leading" velocity profiles produced by rotating disks 
(Fig. 3 in Prochaska & Wolfe 1997). Circle is counter-clockwise rotat- 
ing disk with Vrot = 250 km s^^. Dots are intersection points between 
sightlines and midplane of disk. Vertical dashed line is sightline for in- 
tersection points 1 and 2, and horizontal line is kinematic major axis. 
Absorption profiles produced at indicated locations 



The evidence for rotating disks stems from accurate velocity profiles of 
metal lines that J. X. Prochaska and I have acquired with HIRES, the echelle 
spectrograph on the Keck I telescope. Along with profiles kindly supplied by 
W. L. W. Sargent, we have acquired accurate profiles for low-ions such as Fe"*" 
and Si"*" for over 40 high-z damped Lya systems. Low ions are crucial as they 
trace the neutral gas which gives rise to damped Ijya lines. 

The velocity profiles (see Wolfe 1999; Prochaska & Wolfe 1997) consist of 
multiple narrow-velocity components spanning velocity intervals Av between 20 
and 290 km s^^. The components are not randomly distributed in velocity space. 
Rather, the component with the largest optical depth is at the profile edge in 75 
% of the profiles with At; > 40 km s~^. The most straightforward explanation 
for this "edge leading" asymmetry is that the damped Lya sightlines traverse 
rotating disks in which the density of absorbing clouds decreases exponentially 
in the radial and vertical directions with scale lengths Ra and h (Prochaska & 
Wolfe 1997). Figure 1 shows how "edge- leading" profiles naturally arise from 
randomly selected impact parameters. The figure also shows why Av decreases 
with increasing impact parameter. While other mechanisms have been suggested 
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Figure 2. [Fe/H] versus Av for damped Lya systems at (a) z < 3 
and (b) z > 3, where [Fe/H]=log(Fe/H)-log(Fe/H)0 

to explain the profile asymmetries (e.g. Haehnelt et al. 1997), exponential 
rotating disks provide the most straightforward explanation. 

Using Monte Carlo techniques, we tested the semi-analytic models of Mo et 
al. (1998) in which centrifugally-supported disks are at the centers of dark- 
matter halos with mass distribution determined by Press-Schecter functions 
and CDM power spectra. We generated synthetic velocity profiles by letting 
randomly selected sightlines traverse randomly-oriented exponential disks with 
rotation speeds Vrot equal to yjGM{r)/r at the virial radius. We found the Af 
predicted by single-disk CDM models cannot account for the significant fraction 
of Au between 100 and 300 km s~^, since the median Vrot predicted by single- 
disk CDM models are less than 100 km s~^. However, the TF model in which 
the input distribution of Vrot at high redshift is determined by the Tully-Fisher 
relation and the Schecter luminosity function is compatible with the data owing 
to the larger fraction of disks Vrot > 100 km s~^ (Wolfe & Prochaska 2000). 
Thick disks with /i~ 0.3Rd are required to obtain Av > 100 km s~^. 



3. Metals in Damped Lya Systems 

3.1. Metallicity versus Kinematics in Damped Lya Systems 

We also used the HIRES velocity profiles to determine the metallicities of damped 
Lya systems (Prochaska & Wolfe 2000). Combining metallicity and kinematic 
information we find a phenomenon that bears upon the presence of rotating 
disks at high redshifts. And this is the clear deficit at z < 3 of metal-poor 
systems with Av > 100 km s~^ shown in Figure 2a. The absence of these sys- 
tems is real and not an artifact due to selection effects. This is indicated by 
the presence of such systems in the z > 3 sample (Fig. 2b). Moreover, dust 
would lead to the absence of metal-rich rather than metal-poor systems. The 
enlarged sample of [Fe/H] abundances strengthens the case first made with a 
smaller sample of [Zn/H], Av pairs (Wolfe & Prochaska 1998). We interpret 
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Figure 3. Abundance ratios versus [Fe/H] for damped systems (cir- 
cles) and thick disk stars (stars) (from Prochaska et al. 2000) 



this effect to arise from negative radial gradients in [Fe/H] in rotating disks. 
High [Fe/H] are required at small radii because low impact parameters, b, are a 
necessary condition for large Av (see Fig. 1). We find a wider distribution of 
b contributes to systems with smaller Av. Thus, at small Av, the lower [Fe/H] 
come from low [Fe/H] at large radii, while higher [Fe/H] again come from small 
b but in disks that are more face-on. Therefore, the increase in scatter of [Fe/H] 
with decreasing Av is naturally produced by the rotating disk models. However, 
since variations in metallicity in nearby spirals rarely span more than 1 dex, a 
positive correlation between [Fe/H] and M would help to explain the full effect. 



3.2. Abundance Patterns and the Thick Disk Connection 

Figure 3 plots abundance ratios of various elements versus [Fe/H] for damped 
Lya systems and stars in the Galaxy thick disk. Comparison with the thick disk 
is motivated by the large vertical scale-heights required in our model gaseous 
disks, since this raises the possibility that high-z damped Lya systems are the 
progenitors of thick stellar disks. To demonstrate the plausibility of this hy- 
pothesis, we showed that the metal-rich damped Lya systems contain sufficient 
baryons at mean metallicity, < [Fe/H] > = — 0.6, to account for the the mass 
content and metallicity of thick disks (Wolfe & Prochaska 1998). Here I discuss 
further evidence in support of this idea. 
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The abundances in Figure 3 were inferred from HIRES data obtained for 
a kinematically selected sample of stars in the thick disk and our sample of 
damped Lya systems. The stellar atmospheres programs used to deduce the 
thick disk abundances are described in Prochaska et al. (2000). The figure 
plots results for a , Fe peak, and light elements. As [Fe/H] increases we find 
tentative evidence that [X/Fe] for the damped systems evolves toward [X/Fe] 
for the thick disk when X = Si, Zn, and Mn. The results are especially striking 
for [Mn/Fe]. We find two possible interpretations for these trends. The first is a 
nucleosynthetic explanation. Lu et al. (1996) argued that the damped systems 
exhibit classic type II Sn patterns with enhanced [a/Fe] ratios, deficient [Mn/Fe] 
ratios (due to the odd-even effect), and unpeculiar [Cr/Fe] and [Ni/Fe] ratios. 
The problem with this argument was that [Zn/Fe] = in stars with a wide range 
in metallicities; whereas [Zn/Fe] ~ 0.3 in damped systems. With our more 
accurate spectroscopy we find [Zn/Fe] = 0.1 in thick-disk stars where [Fe/H] 
= — 0.6. And there is tentative evidence that [Zn/Fe] is higher for thick-disk 
stars with lower [Fe/H] . That Zn does not track Fe is not surprising since recent 
calculations indicate Zn is not produced in the same type I Sn responsible for 
most of the Fe. Rather, Zn may be produced in neutrino-driven winds in type 
II Sn (Hoffman 1996)). Therefore, it is likely that the [Zn/Fe] ratio in damped 
Lya systems is partially nucleosynthetic in origin. On the other hand, Pettini 
et al. (1997) argued that the abundance patterns in damped systems are caused 
by depletion of heavy elements on to dust grains. This is consistent with all 
the observed patterns except for [Mn/Fe] and possibly [Ti/Fe]. The strongest 
argument in favor of this explanation was the high [Zn/Fe] ratio. But we have 
shown that part of this ratio probably originates from nucleosynthesis, and so a 
more likely explanation is that the damped Lya abundance pattern results from 
small amount of dust superposed on a type II Sn abundance pattern. 

As a result, we conjecture that metal-poor damped Lya systems evolve 
into metal-rich objects that become thick stellar disks at high redshifts, say z 
~ 3. Since the velocity dispersion of the model gaseous disks, a ~ 10 km s^^, 
subsequent merger events are required to heat the stars to the observed velocity 
dispersion of the thick disk, cr ~ 40 km s^^. In that sense our model does not 
differ substantially from the scenario put forward by Wyse (these proceedings). 



4. Connecting Lyman Break Galaxies to Damped Lya Systems 

For the past two years, our group (Gawiser, Prochaska, Cooke, & Wolfe) has been 
searching for Lyman break galaxies in fields centered on quasars with foreground 
damped Lya systems. Our goals are to (i) find galaxies physically associated 
with the damped systems, and (ii) determine how often the damped Lya red- 
shifts coincide with the redshift spikes characterizing the distribution of known 
Lyman break galaxies (Steidel 1999). The correlation between damped systems 
and Lyman break galaxies is very interesting because it tells us about the bias 
of the damped Lya galaxies, and, as a result, about their mass. Stated differ- 
ently, if damped Lya redshifts are preferentially located in Lyman break redshift 
spikes, this would mean the amplitudes of density fluctuations, {6p/p)DLA, des- 
tined to become damped Lya galaxies are too low to collapse and form galaxies 
unless they are boosted by the enhanced mass density provided by the Lyman 
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break clusters. This would indicate higher-than-typical masses given the nature 
of the power spectrum. 

Our survey initially focused on B drop outs as there is no U band photom- 
etry currently available on the Keck telescopes. This restricted the sample to 
damped Lya systems with z =(3.8,4.5). We have acquired B,R,I images for 
about four damped Lya fields in which the la sensitivies correspond to to -B ~ 
29.5, R ^ 26.5, and / « 26.0. As a result we are able to detect blue drop-outs 
down to i? ~ 26.0. We select candidate blue drop-outs using criteria established 
by Steidel and then acquire multi-object spectra for about 20 candidates per 
setting using the LRIS multi-object spectrograph. The data have been reduced 
for one field containing damped systems at z = 3.871 and 4.072. While there 
are no confirmed galaxies near the high redshift systems, we have identified four 
galaxies within Az = 0.02 of the lower redshift system. That is, four out of 
a total of twelve identified galaxies lie within a Steidel-like redshift spike that 
includes the damped redshift. If future observations show this this effect to be 
generic, then some fraction of damped systems may evolve into galaxies likely 
to be found in cluster environments, such as ellipticals. 



5. [C II] 158 nm Emission from Damped Lya Systems 

[C II] 158 fira emission results from transitions between the P3/2 and P^/2 fine 

structure states in C'^. It dominates cooling by the Galaxy ISM with a lumi- 
nosity L([C II]) =5xlO'^L0 (Wright et al. 1991). Significantly, most of the 158 
fiiQ emission from the ISM in the Galaxy and in other spiral galaxies where it 
is detected comes from the cold neutral medium (CNM) rather than the warm 
neutral medium, star- forming regions in spiral arms (Madden et al. 1993), or 
PDRs on the surfaces of molecular clouds. The last point is especially relevant 
for damped Lya systems where molecules are rarely detected (cf. Lu et al. 
1999). 

The 158 ^m luminosities of damped Lya systems can be estimated from 
A^(C II*), the column density of C+ ions in the P3/2 state, which is inferred 

from the absorption profiles of the UV transition C II* 1335.708. The 158 /im 
luminosity per H atom is given by 

Ic = /iz/2iiV(Cir)A2i/7V(HI)ergs-iH"\ (1) 

where hi/21 is the energy of the transition and A21 is the coefficient for spon- 
taneous photon emission (Lu et al. 1999). Since -^([C II])=</c>Af (H l)/mH, 
where <lc> is the density-weighted average of Ic, we can estimate the C II 
luminosity from damped systems provided we know the H I mass and <lc>- 

In principle we can determine <lc> from the distribution of Ic measured 
along individual sightlines. Figure 4 shows estimates of Ic for 11 damped Lya 
systems. We also show Ic deduced for separate sightlines through the ISM and 
<lc> estimated for the ISM. We compute <lc> for the damped sample by as- 
suming exponential disks. In that case <lc> ~ Ic for impact parameters b ~ 
2R(i. Our Monte Carlo simulations show that most damped sightlines have b < 
(4— 5) xi?(i. As a result <lc> corresponds to largest (2/4.5)^ x 11 «:: 2—3 values of 
Ic in Figure 4. We estimate log<^c> = —26.6 ±0.2. While somewhat uncertain. 
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Figure 4. [CII] 158 /im emission per H atom versus N(HI). Circles 
are damped Lya systems with upper limits and one lower limit shown. 
Small stars are estimates for separate ISM sightlines. Large star is 
global average of ISM from COBE 



<lc> for the damped systems is clearly more than 30 times lower than for the 
ISM. Since the [C II] 158 /im line is the dominant coolant for the CNM, <lc> 
equals the heating rate for steady state conditions. As a result, the low value 
of <lc> indicated for damped Lya systems suggests low heat inputs into the 
protogalactic gas at high redshifts. 

We also use <lc> to compute the distribution of L([C II]) of damped Lya 
systems in various cosmogonies. Mo et al. (1998) predict the fraction of inter- 
cepted damped Lya systems with dark-halo masses exeeding M. We convert 
M to 158 fim emission noting that L{[C Il]) = <lc>m^M / firriH , where rrid is the 
fraction of dark halo mass in the disk. Assuming log</c>=— 26.5 erg s~^II~^, /i 
= 1.4, and md=0.05, we find L([C II])=2xW'^ {M/W^'^Mq)Lq. We compute the 
flux density at a given redshift by assuming a rectangular velocity profile with 
width = 2xV^o(sin(z), where sin(i) = 0.5. The results for the ACDM and TF 
cosmogonies are shown in Figure 5. At z ~ 3 the predicted ALMA detection rate 
is < 3 % the CDM model, and about 25 % in the TF model. Thus detection 
of [C II] 158 /im emission in high-z damped systems would be highly signifi- 
cant. Moreover, measurements of the velocity profile of the emitting gas could 
distinguish rotating disks from protogalactic clumps (Haehnelt et al. 1997). 
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Figure 5. Fraction of damped Lya systems detected with 158 fim 
flux density exceeding Si,, (a) Curves are predictions for ACDM model 
at redshift z = 1,2,3,4. Error bars are 3o" ALMA sensitivities for S,y at 
each redshift (Brown 2000). (b) Same as (a) except for TF model 
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